function [logL]=likelihood_obj(y,R)

alpha=y(1); 
alpha_2=y(2);
lambda_4=y(3); 
lambda_5=y(4); 
kappa =y(5);   
ln_sigma_v =y(6); 

L = integral(@(v) likelihoods(y,v,R).*normpdf(v,0,exp(ln_sigma_v)),-Inf,Inf,'ArrayValued',true);

logL = -sum(log(L));
